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HIGH  ORDER  "ZIP”  DIFFERENCING  OF  CONVECTIVE  TERMS 

I.  INTRODUCTION 

Consider  equations  of  the  form 

wt  +  fx  -  0  (1) 

where 

f  =  w  +  f  (2) 

and  hence 

w  +  (wv)  +  f*  *  0  (3) 

t  x  x 

Here  w,  f  and  v  are  functions  of  the  independent  variables  x  and  t. 

The  second  term  in  Eq.  (3)  is  called  a  convective  term  and  is  the 
subject  of  this  paper.  We  shall  say  that  a  finite  difference  approxi¬ 
mation  to  Eq.  (3)  is  in  conservation  or  "flux"  form  when  the  approxi¬ 
mation  can  be  written 

Wi+1  *  wi  -  ^  ^i"1  tFi+Ss  -  (4) 

Here  the  subscripts  refer  to  the  spatial  grid  points  x^,  the  super¬ 
scripts  to  the  temporal  grid  points  tn,  At  m  tn+^  -  tn,  and  Ax^  s 

T  (x1+1  -  Xj_^).  We  shall  assume  heneceforth  that  Ax^  is  independent 
of  i,  and  denote  the  quantity  simply  by  Ax.  In  this  author's  opinion, 
nonuniform  meshes  are  best  handled  by  coordinate  transformation  or 
mapping  techniques[3,4,13].  The  ^  are  called  transportive  fluxes, 
and  are  functions  of  the  f^  at  one  or  more  of  the  time  levels  tn.  For 
a  given  time  level,  the  functional  dependence  of  F  on  f  can  be  written 
to  achieve  any  desired  order  of  spatial  accuracy  (see  £l.  Appendix]). 
Manuscript  submitted  March  3, 1980 
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For  instance,  centered  second  order  differences  require 


Fi4** 


+  fi> 


(5) 


while  centered  fourth  order  differences  give 


Fi-Hs  "  12  (fi+l  +  V  “  12  (fi+2  +  fi-l) 


(6) 


For  the  purposes  of  this  paper  we  wish  to  divide  the  net  transportive 
flux  F^ p  into  two  components:  |^,  the  convective  component,  and 

F^(^,  the  non-convective  component: 


Fi+*S  =  Fi4is  + 


(7) 


Thus  F®^  and  F^^  are  the  fluxes  corresponding  to  the  second  and 
third  terms  respectively  of  Eq.  (3).  We  shall  consider  only  F^^ 
here. 

II.  SECOND  ORDER  ZIP  FLUXES 

A  straightforward  implementation  of  Eq.  (5)  for  the  convective 
flux  would  yield 


Fi+>S  “  2  ^Wi+1  Vi+1  +  Wi  Vi^  (< 

However,  as  is  pointed  out  in  the  classic  1968  paper  by  Hirt  [2],  a 
better  form  for  this  convective  flux  is 


(9) 


14*5 


1  Cvi+1  Wi  +  Wi+1  Vi] 


Hirt  refers  to  this  form  as  "ZIP"  differencing,  a  convention  we  shall 
keep  here.  Hirt's  "heuristic"  stability  analysis  of  difference  equations 
entails  expanding  each  term  in  the  difference  equation  in  a  Taylor 
series  to  obtain  the  differential  equation  one  is  actually  solving. 
Examination  of  the  properties  of  this  new  differential  equation  sheds 
much  light  on  the  stability  and  error  characteristics  of  the  difference 
equation.  Expanding  our  flux  form  representation  for  the  second  term 
in  (3)  we  find 


(tr) FD  *  -1  ^  +  « 


(10) 


where  the  superscript  refers  to  the  finite  difference  approximation 
and  TE  represents  the  truncation  error  terms. 

For  Eq.  (8)  we  find  that 


TE  -  |  (  Ax2  +  )  Ax4  +  0  (Ax6) 

dx  i  dx  i 


while  for  the  ZIP  flux  (9)  we  find 


(ID 


1 


In  Hirt's  heuristic  stability  analysis  of  (3),  one  looks  for  truncation 

b2 

error  terms  of  the  form  — ^  since  such  terms  will  destabilize  numerical 
solutions  to  Eq.  (1)  if  they  are  negative  and  artificially  damp  the 
solutions  if  they  are  positive.  It  should  be  noted  that  for  certain 
explicit  non-dissipative  schemes  (leapfrog,  for  example)  terms  of  this 
form  are  destabilizing  regardless  of  their  sign.  A  look  at  (12)  shows 
that  terms  of  this  fotm  are  simply  not  present.  One  can  extend  Hirt's 
argument  and  look  for  terms  involving  any  even  derivatives  of  w  with 
respect  to  x,  since  any  such  terms  could  potentially  be  destabilizing, 
but  these  too  are  absent  when  ZIP  fluxes  are  used.  The  odd 
order  derivatives  in  the  truncation  error  terms  will  give  rise  to 
error,  but  these  will  be  of  a  dispersive  rather  than  of  a  dissipative 
or  antidissipatlve  nature. 

By  contrast,  if  we  examine  the  truncation  errors  given  by  (11), 
we  find  that  the  leading  term  is 


1 

6 


r  i.  <^w  dv  .  1  2&v  .  1  djy  _1 

2  dx2  *  2  dx2  *  6  dx3  6 


(13) 


which  will  contribute  a  destabilizing  or  dissipative  term  except  in  the 
dv 

trivial  case  *  0.  Thus  although  both  forms  for  the  convective  flux 
have  the  same  formal  order  of  accuracy  in  terms  of  Ax,  we  find  the  ZIP 
form  preferable  from  stability  considerations. 
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III.  HIGH  ORDER  "ZIP"  FLUXES 


The  extension  of  the  ZIP  concept  to  higher  order  can  be  seen  by 
FD 

expanding  (wv)x  for  the  second  order  ZIP  flux  (9) 


FD 

^  -rCi-,3 


1 3x 


i  [»1+1  v1  +  »±  »1+1  - »  vr^-i  v±3 


i  i*"1  »i+i  -  ”i-i]  vi  +  T  i*"1  [v1+i  -  vi-i  ]  wi 


VFD 


\FT> 


+  w 


1  i^J 


a*) 


So  we  see  that  ZIP  differencing  of  the  convective  flux  is  equivalent  to 

taking  a  central  second  order  finite  difference  approximation  to  the 

term  wv  +  vw  rather  than  (wv)  .  We  point  out  that  Cheng  and  Shubin 
XX  X 

[5]  have  also  noted  that  second  order  central  differencing  of  the  term 
wvx  +  vwx  leads  to  a  conservative  scheme,  something  we  have  proven  here 
by  showing  the  ZIP  form  to  be  the  equivalent  convective  flux.  It 
would  seem,  then,  that  higher  order  forms  for  ZIP  fluxes,  if  they 
exist  at  all,  are  equivalent  to  higher  order  central  difference 
approximations  to  the  form  wvx  +  vwx.  In  fact,  it  is  now  quite  clear 
why  the  ZIP  fluxes  work  as  they  do:  For  central  differences  of  all 

orders  of  accuracy,  the  truncation  errors  of  approximations  to 
first  derivatives  involve  only  odd  order  derivatives  of  the  function. 
Hence,  the  form  wvx  +  vwx  can  never  produce  terms  involving  even  order 
derivatives  of  w,  as  long  as  centered  finite  differences  are  used. 
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By  way  of  contrast  "straightforward"  central  differencing  of  the  form 

(wv)  will  still  produce  only  odd  order  derivatives,  but  those 
x 

derivatives  will  operate  on  the  quantity  wv,  and  can  always  be  ex¬ 
panded  to  show  that  terms  involving  the  undesirable  even  order  de¬ 
rivatives  of  w  are  present.  For  example,  the  "straightforward"  fourth 
order  convective  flux 

Fl*  •  it  (Vi  +  Vl  vi+l>  •  it  <U1-1  Vi-1  +  "i+2  vl+2>  (15> 


gives  the  truncation  error 


TE  -  -  gj  Ax4  +  0  (Ax6) 

3x 


(16) 


Expansion  of  the  leading  term  shows  components  involving  both 
a2  zfi 

ox  dx 

Without  further  ado,  we  give  here  the  appropriate  ZIP  fluxes 
for  the  indicated  spatial  order  of  accuracy: 

Second  Order: 


i+*S 


T  l>  v  +  w.  v  ] 


i+1  i 


i  i+1 


(17) 


Fourth  Order: 

rl+>s  -ft  wi+l  V1  +  “l  Vl+1] 

-  V1  +  *1  vl+2  +  wl+l  vl-l  +  ”1-1  W  (1») 
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Sixth  Order: 


Frt  ‘  I  Csi+i  wi  +  ”i  vi+i^ 


20 


[wi+2  v±  +  wt  vi+2  +  wi+1  vi_1  +  w1_1  v±+13  (19) 


+  6^  Cwi+3  vi  +  Wi  Vi+3  +  wi+2  Vi-1  +  Wi-1  Vi+2 


+  Wi+1  Vi-2  +  Wi-2  vi+l] 


In  general 


c  _  J  1  rN 
i+*5  k  k  ki 


(20) 


where  N  =  order  of  approximation  (even)  -5-  2 

_N  N 
1  "  NH 


rN  -  pN  N  -  k  +  1 
k  '  k-1  *  N  +  k 

k-1 

Pki  =  (wi-j+k  Vi-j  +  Vj  Vi-j+k} 


(21) 

(22) 


(23) 


The  reader  can  verify  easily  that  the  quantity  [F^  ^  -  F^_^]  Ax  ^  does 
in  fact  reproduce  the  appropriate  order  finite  difference  approximation 
to  the  quantity  (w  v^  +  wx  v)^  when  the  above  expressions  are  used. 

TV.  NUMERICAL  EXAMPLE 

We  choose  as  our  computational  example  the  exploding  diaphragm  or 
Riemann  problem  for  the  one-dimensional  equations  of  ideal  inviscid 
compressible  fluid  flow: 
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p 

pv 

pE 


0 


L 


t 


+ 


*  Pv 

2 

pv 

pEv 


x 


+ 


P 

Pv 


=  0 


X 


(24) 

where  p,  v,  P  and  E  are  the  fluid  density,  velocity,  pressure  and 
specific  total  energy  respectively.  Also, 

P  3  (y-1)  (pE  =  \  pv^)  (25) 

where  we  will  choose  Y=  1.4.  Our  initial  conditon  consists  of  two 
constant  states,  one  to  the  left  and  one  to  the  right,  separated 
by  a  discontinuity  which  is  assumed  to  lie  midway  between  two  grid 
points,  in  this  case  points  50  and  51.  For  these  calculations  * 

P,  c.  =  1.0,  P  ,  =  10.0,  and  P  .  =  100.0.  All  calculations  were 

done  on  a  uniform  mesh  of  100  grid  points  using  a  Courant  number  of 
0,1,  and  utilizing  the  second  and  fourth  order  leapfrog-trapezoidal 
flux-corrected  transport  (FCT)  algorithms  described  in  [l].  For  sim¬ 
plicity  the  Boris-Book  flux  limiter  [l,6j  was  used.  Briefly,  FCT 
algorithms  solve  Eq.  (1)  by  forming  a  point  by  point  weighted  average 
of  two  transportive  fluxes,  one  chosen  to  give  monotone  results  (free 
of  spurious  oscillations)  for  the  problem  at  hand,  and  the  other  to 
give  formally  high  order  accuracy.  The  high  order  fluxes  are  weighted 
as  heavily  as  possible  without  violating  monotonicity  constraints,  a 
process  referred  to  as  "flux  limiting."  In  this  example  our  monotone 
flux  is  given  by  the  first  order  Rusanov  scheme  [7],  while  our  high 
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order  fluxes  are  given  by  the  aforementioned  second  and  fourth  order 
leapfrog  trapezoidal  schemes.  The  tests  we  will  make  of  ZIP  versus 
non-ZIP  differencing  of  convective  terms  refer  only  to  the  high  order 
portion  of  the  FCT  algorithm.  FCT  algorithms  are  extremely  effective 
in  suppressing  the  numerical  oscillations  and  instabilities  which 
would  otherwise  occur  with  the  high  order  scheme.  Because  of  this, 
there  may  be  a  tendency  on  the  part  of  the  user  to  believe  that  great 
care  is  not  required  in  choosing  the  higher  order  portion  of  the  total 
FCT  scheme  since  the  flux  limiting  will  "save"  him.  We  shall  show  that 
this  is  not  always  the  case. 

We  have  written  Eq.  (24)  to  show  explicitly  the  convective  terms 
as  the  second  terms  in  the  (vector)  equation.  In  Fig.  1  we  show  the 
calculated  (data  points)  and  analytic  (solid  line)  density  profiles 
after  500  time  steps  using  the  second  order  accurate  convective  fluxes 
given  by  the  "straightforward"  evaluation  Eq.  (8)  and  by  the  ZIP 
evaluation  Eq.  (17).  We  note  that  analytically  we  have  a  shock  wave 
moving  to  the  left,  followed  by  a  contact  discontinuity,  and  a  rare¬ 
faction  fan  moving  to  the  right.  The  straightforward  calculation  is 
seen  to  be  afflicted  by  several  large  scale  numerical  oscillations  while 
the  ZIP  calculation  yields  a  reasonably  accurate  and  oscillation-free 
profile,  even  though  the  implementation  of  FCT  has  kept  both  calculations 
stable.  From  (3),  (10),  and  (13)  we  see  that  straightforward  differ¬ 
encing  is  equivalent  to  adding  a  diffusion  term  with  coefficient 
-(dv/dx)/2  to  the  right  hand  side  of  (3).  When  dv/dx  is  positive,  as 
in  the  rarefaction  fan,  we  expect  numerical  instability;  while  in  the 
shock  (dv/dx  <  0)  we  expect  artificial  dissipation.  We  note  that  all 


9 


of  the  regions  of  roughness  or  oscillation  in  Fig.  la  represent  fluid 
elements  that  are,  or  at  one  time  were,  in  the  rarefaction  fan;  and 
further  that  the  density  jump  at  the  shock  is  smaller  and  occupies  a 
greater  number  of  cells  than  in  the  corresponding  ZIP  calculation. 

Fig.  lb.  The  ZIP  calculation  is  obviously  far  superior.  To  be  fair, 
other  procedures  could  have  "fixed"  the  problem.  In  fact  our  choice 

of  the  virtually  non-dissipative  leapfrog-trapezoidal  scheme  is 
probably  not  a  wise  one  for  a  calculation  of  this  type,  where  discon¬ 
tinuities  abound.  Nonetheless  we  do  feel  this  example  makes  a  valid 
comparison  of  convective  differencing  schemes. 

In  Fig.  2  we  show  the  same  comparison  as  in  Fig.  1,  but  this  time 
for  the  fourth  order  "straightforward"  flux  Eq.  (15)  and  for  the  fourth 
order  ZIP  flux  Eq.  (18).  Again  the  "straightforward"  calculation 
gives  numerical  oscillations  whereas  ZIP  differencing  experiences  no 
apparent  difficulties.  Comparing  Figures  1  and  2  we  see  an  improvement 
in  the  solution  as  we  go  from  second  to  fourth  order,  a  behavior 
one  would  expect  based  on  a  linear  dispersion  analysis.  Further 
improvements  in  the  ZIP  solutions  are  seen  as  the  order  of  the 
approximation  is  increased  to  sixth  and  higher,  but  the  oscillations 
associated  with  the  "straightforward"  differencing  remain  intact. 

In  fact  the  improvement  shown  here  is  quite  mild  compared  to 
what  we  have  seen  in  other  types  of  calculations  as  we  increase 
the  spatial  order  of  the  approximation.  It  now  appears  that  the 
"ultimate"  FCT  scheme  will  be  either  a  very  high  order  polynomial-based 
approximation  scheme  or  a  pseudospectral  scheme  (see  [8]).  Both  the 
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second  and  fourth  order  calculation  smear  the  contact  discontinuity 
more  than  is  desirable.  This  is  due  to  the  overly  conservative 
assumptions  of  the  flux  limiter  used  (the  "clipping"  problem  -  see  [l]) 
at  early  times,  and  can  be  remedied  by  using  a  different  flux  limiter 
[l]  and  using  monotonicity  constraints  whose  description  is  beyond  the 
scope  of  this  paper. 

Before  leaving  this  section  we  briefly  present  several  further 
calculations  which  may  be  of  some  interest  to  the  reader. 

Another  very  popular  way  of  writing  a  second  order  convective 


flux  is 


+ 


(vi+l 


V 


(26) 


In  fact  MacCormack  [14]  has  also  concluded  that  the  straightforward 
flux  (8)  is  to  be  avoided  in  certain  regions  and  is  to  be  replaced  by 
(26)  there.  Comparison  of  (8),  (9),  and  (26)  shows  that  (26)  is 
nothing  more  than  an  average  of  second  order  ZIP  and  straightforward 
fluxes.  The  results  for  our  Riemann  test  problem  using  (26)  for  the 
convective  fluxes  are  shown  in  Fig.  3.  We  note  that  the  results  are 
far  superior  to  those  for  the  straightforward  fluxes.  Fig.  la,  but 
not  quite  as  good  as  those  for  the  ZIP  fluxes.  Fig.  lb. 

We  implied  earlier  that  some  of  the  problems  associated  with 
straightforward  fluxes  could  be  overcome  by  using  a  dissipative  scheme. 
In  Fig.  4,  we  show  the  results  of  using  the  fourth  order  straightforward 
flux  (15)  plus  a  dissipative  flux  of  the  form 
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-  3T  <  I 


vi+i  +  vi, 


ci+l  +  ci 


)  [  w 


i+2 


'  “l-l  -3«»i+l 


-  wt)  ] 


(27) 

which  approximates  a  fourth  order  dissipation  term  involving  the  fourth 
derivative  of  v.  Here  c^  is  the  sound  speed  at  grid  point  i.  Note 
that  the  results  are  much  improved  over  Fig  2a.  In  Fig  5  we  show  the 
results  of  adding  the  same  dissipative  flux  (27) to  the  fourth  order 
ZIP  flux  (18).  The  results  are  degraded  only  slightly.  The  idea  of 
using  fourth  order  dissipation  with  non-dissipative  fourth  order  schemes 
is  originally  due  to  Kreiss  and  Oliger  [15],  and  the  obvious  extension 
to  even  higher  order  schemes  is  the  subject  of  study  by  the  present 
author. 

The  reader  may  note  that  in  the  energy  component  of  (24),  the 
Pv  term  could  easily  have  been  treated  using  the  ZIP  format.  In  fact, 
by  (25)  P  is  proportional  to  one  component  of  pE  and  is  therefore  sub¬ 
ject  to  the  same  kind  of  nonlinear  instabilities  as  are  the  other  con¬ 
vective  terms.  In  Fig  6  we  show  the  results  of  using  the  fourth  order 
ZIP  flux  (18)  for  both  the  pEv  and  the  Pv  terms  in  the  energy  equation. 
Although  there  is  very  little  difference  in  this  case  from  the  straight¬ 
forward  treatment  of  Pv,  the  ZIP  flux  form  is  recommended  on  theo¬ 
retical  grounds. 

V.  CONCLUSIONS 

We  have  shown  that  ZIP  differencing  of  the  term  (wv)^  is  the 
equivalent  of  a  centered  differencing  of  the  alternative  form 
wv^  +  w^v,  both  in  the  second  order  limit  and  with  regard  to  the 
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kinds  of  truncation  error  terms  that  are  produced  for  the  higher 
orders.  Explicit  expressions  for  the  calculation  of  arbitrarily  high 
order  ZIP  fluxes  have  been  given,  and  computational  examples  have 
been  given  to  show  the  advantages  of  the  ZIP  flux  form  over  the 
"straightforward"  flux  form  within  the  context  of  FCT  algorithms. 

It  seems  to  us  that  the  "heuristic"  stability  analysis  of  Hirt  [2] 
has  again  proved  to  be  a  reliable  tool  in  analyzing  numerical  schemes. 
We  should  point  out  to  the  reader  that  the  ZIP  flux  concept  is  but 
one  example  of  a  class  of  numerical  techniques  based  on  Hirt's  analysis 
known  as  truncation  error  cancellation  (TEC)  which  are  now  quite  highly 
developed  and  which  have  even  been  automated  [9-11].  Similar,  but 
not  identical,  techniques  have  been  given  by  Warming  and  Hyett  [12], 
Lerat  and  Peyret  [16],  Yanenko  and  Shokin  [17],  and  by  Cheng  [18], 
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2ND  DRDER  STRAIGHT 


X 

Fig.  la  -  Comparison  of  analytic  and  computed  density  profiles  (solid 
lines  and  data  points  respectively)  for  the  exploding  diaphragm  problem 
given  in  the  text  for  second  order  straightforward  convective  fluxes 
Eq.  (8) 
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Pig.  lb  -  Comparison  of  analytic  and  computed  danaity  profiles  (solid 
lloes  and  data  points  respectively)  for  the  exploding  diaphragm  problem 
given  in  the  text  for  second  order  ZIP  fluxes  Eq.  (17) 
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4TH  DRDER  STRAIGHT 
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Fig.  2a  -  Same  as  Fig.  1  but  for  fourth  order  straightforward 
convective  fluxes  Eq.  (15) 
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4TH  BTR  D-l/12 


Fig.  4  -  Same  as  Fig.  2a  but  with  the  addition  of  a  fourth  order 
dissipation  term  given  by  Eq.  (27) 
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DENSITY 
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5  -  Same  as  Fig.  2b  but  with  the  addition  of  a  fourth  order 
dissipation  term  given  by  Eq.  (27) 
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01CY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

SAM)  I A  LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  Ml  87115 

01CY  ATTN  OOC  CON  FOR  J.  MARTIN 

01CY  ATTN  DOC  CON  FOR  W.  BROWN 

01CY  ATTN  DOC  CON  FOR  A.  TH0RN8R0UGH 

01CY  ATTN  DOC  CON  FOR  T.  WRIGHT 

01CY  ATTN  OOC  CON  FOR  D.  DAHLGREN 

01CY  ATTN  DOC  CON  FOR  5141 

01CY  ATTN  DOC  CON  FOR  SPACE  PROJECT  OIV 


SAM)I A  LABORATORIES 
LIVERMORE  LABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  B.  (RJRPMEY 
01CY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  D.  GALE 


OTHER  GOVERMCNT 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI,  RM  5G48,  HQ  BLDG 
WASHINGTON,  D.C.  20505 

01CT  ATTN  OSI/PSIO  RM  5F  19 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STAMOARDS 
WASHINGTON,  O.C.  20254 

CALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
01CY  ATTN  R.  MOORE 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMfONICATlONS  t  INFO  ADMIN 
BOULDER,  CO  80505 

01CY  ATTN  A.  JEAN  (UNCLASS  0N.Y) 
01CY  ATTN  W.  UTLAUT 
01CY  ATTN  0.  CROMEIE 
01CY  ATTN  L.  BERRY 


NATIONAL  OCEANIC  8  ATMOSPHERIC  AOMIN 
ENVIROMPCNTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80502 

01CY  ATTN  R.  GRUBB 

01CY  ATTN  AERONOMY  LAB  6.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.  0.  BOX  92957 
LOS  ANGELES,  CA  90009 

01CY  ATTN  I.  GARFUMCEL 
OICT  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  N.  STOCKWELL 
01CY  ATTN  D.  OLSEN 
01CY  ATTN  J.  CARTER 
01CY  ATTN  F.  MORSE 
01CY  ATTN  SMFA  FOR  PWW 

ANALYTICAL  SYSTEMS  ENGINEERING  CORF 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01805 

01CY  ATTN  RADIO  f:iB«ES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.  0.  BOX  985 
BERKELEY,  CA  94701 

01CY  ATTN  J.  WORKMAN 

BOEING  COMPANY,  THE 
P.  O.  BOX  5707 
SEATTLE,  WA  98124 

01CY  ATTN  G.  KEISTER 
OVC!  ATTN  0.  (ARRAY 
01CY  ATTN  G.  MALL 
01CY  ATTN  J.  KEMNEY 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
IPAPS,  B-019 
LA  JOLLA,  CA  92095 

01CY  ATTN  HEMtY  G.  BOOKER 

BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  55807 

01CY  ATTN  ROMEO  A.  OELIBERIS 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02159 

01CY  ATTN  0.  B.  COX 
01CY  ATTN  J.  P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  BLVO 
FALLS  CHURCH,  VA  22046 
01CY  AT^N  M.  BLANC 
01CY  ATTN  JOHN  SPOOR 
01CY  ATTN  C.  NAIL 


COMSAT  LABORATORIES 
LINTH1CUM  ROAD 
CLARKSBURG,  MO  20754 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGIMCERING 

ITHACA,  NY  14850 

01CY  ATTN  D.  T.  FARLEY  JR 


ELECT *05* ACE  SYSTEMS,  INC. 

BOX  1 359 

RICHARDSON,  TX  750*0 

OICY  ATTN  H.  LOCSTON 

01CT  ATTN  SECURITY  (PAU.  PHILLIPS) 

ESL  INC. 

095  JAVA  MIVE 
SUPMTVALE,  CA  940*6 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAWS  MARSHALL 
OICY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  t  COMMUNICATIONS  CORP 
5939  FABIAN  HAY 
PALO  ALTO,  CA  94)05 

OICY  ATTN  J.  T.  MATT  INGLE Y 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  1)201 

OICY  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  <P.O.  DRAPER  QQ) 

SANTA  BARBARA,  CA  9)102 
OICY  ATTN  DAS  I  AC 
OICY  ATTN  DON  CHANDLER 
OICY  ATTN  TOM  BARRETT 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAFBILL 
OICY  ATTN  MACK  STANTON 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

COURT  STREET 
SYRACUSE,  NY  1)201 

OICY  ATTN  G.  HILLMAN 


GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  9)111 

OICY  ATTN  JOHN  ISE  JR 
OICY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBAMKS,  AK  99701 

CALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CT  ATTN  T.  N.  DAVIS  CUNCL  ONLY) 

OICY  ATTN  NEAL  BROWN  CUNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  OIV 
77  A  STREET 
MCEDHAM,  MA  02194 

OICY  ATTN  MMRSMAL  CROSS 


ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  6180) 

OICY  ATTN  K.  TEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WIGHT  STREET 
URBANA,  IL  60680 

CALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  AMW.YSES 
400  ARKT-NAVY  DRIVE 
ARLINGTON,  VA  22202 

OICY  ATTN  J.  M.  AEIN 
OICY  ATTN  ERMCST  BAUER 
OICY  ATTN  HANS  WOLF HARD 
OICY  ATTN  JOEL  BENGSTON 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MM  01730 

OICY  ATTN  DONALD  HANSEN 

INTO  TEL  t  TELEGRAPH  CORPORATION 
SOO  WASHINGTON  AVENUE 
NUTLET,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMINO  DEL  MAR 
DEL  MAR,  CA  92014 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  TMOMMS  POTEMRA 
OICY  ATTN  JOHV  DASSOULAS 

LOCKHEED  MISSILES  t  SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  W.  L.  IMWP  DEPT  52-12 

KAM5AN  SCIENCES  CORP 
P.  0.  BOX  7465 
COLORADO  SPRINGS,  CO  809)5 
OICY  ATTN  T.  MEAGHER 

LIMKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

M.l.T.  LINCOLN  LABORATORY 
P.  0.  BOX  73 
LEXINGTON,  »M  02175 

OICY  ATTN  OAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGH. IN 
OICY  ATTN  0.  CLARK 

MARTIN  MARIETTA  CORP 
ORLAMCX)  DIVISION 
P.  0.  BOX  5837 
ORLAMCO,  FL  32803 

OICY  ATTN  R.  HEFFNER 


MCOOWCU  00U6US  corporation 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
OICY  ATTN  d.  HOULE 
OICY  ATTN  GEORGE  HROZ 
OICT  ATTN  W.  OLSON 
OICT  ATTN  R.  W.  HALPRIN 
OICT  ATTN  TECWICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
OICY  ATTN  W.  F.  CREVIER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  0.  SAPPEMMELD 
OICY  ATTN  R.  80GUSCH 
OICY  ATTN  R.  HETORICK 
OICY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOULE 
OICY  ATTN  F.  FAdEN 
OICY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGMIRE 
OICY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  0.  BOX  206 
BEDFORD,  MA  01730 

OICY  ATTN  JOM  MORGANS  TERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

westgate  research  park 
1*20  DOLLY  MADISON  BLVO 
MCLEAN,  VA  22101 

OICY  ATTN  W.  MALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVE RFIELO  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E.  C.  FIELD  JR 

PE7P6YL VANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
31*  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16*02 

(NO  CLASSIFIED  TO  THIS  /MJORESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOKTRICS,  INC. 

442  MAR  RETT  ROAD 
LEXINGTON,  MA  02173 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  WA  9*009 

OICY  ATTN  E.  d.  FRfMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

oicy  attn  a.  Thompson 


R  t  D  ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  d.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 


RAM)  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
U  JOLLA,  CA  9203* 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  D.  SACHS 
OICY  ATTN  E.  A.  STRNCER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 
WJNTSVILLE  DIVISION 
2109  W.  a  INTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  J5805 

OICY  ATTN  DALE  H.  01V1S 

SCIENCE  APPLICATIONS,  INCORPORATED 
*400  «STPARK  DRIVE 
MCLEAN,  VA  22101 

OICY  ATTN  d.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
OICY  ATTN  SZ 


SRI  INTERNATIONAL 
533  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  ICILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  L.  L.  COBB 
OICY  ATYN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CMESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEAOABRAM) 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  d.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 


TEOfJOLOCY  INTERNATIONAL  COW 
75  WIGGINS  AVENUE 
8E0FORD,  MA  01730 

01CY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  t  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

RHXROO  BEACH,  CA  90270 

01CY  ATTN  R,  K.  PLE8UCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VISIDY1*,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  IWUSTRIAL  PARK 

BURLINGTON,  HA  01003 

01CY  ATTN  CHARLES  HUMPHREY 
01CY  ATTN  J.  W.  CARPENTER 


